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Abstract 

The BABAR Silicon Vertex Tracker (SVT ) is a five-layer double-sided silicon detector designed to provide precise 
measurements of the position and direction of primary tracks, and to fully reconstruct low- momentum tracks produced 
in e"'"e~ collisions at the PEP-II asymmetric collider at Stanford Linear Accelerator Center. This paper describes 
the design, implementation, performance, and validation of the local alignment procedure used to determine the 
relative positions and orientations of the 340 SVT wafers. This procedure uses a tuned mix of in-situ experimental 
data and complementary lab-bench measurements to control systematic distortions. Wafer positions and orientations 
are determined by minimizing a computed using these data for each wafer individually, iterating to account for 
between-wafer correlations. A correction for aplanar distortions of the silicon wafers is measured and applied. The net 
effect of residual mis-alignments on relevant physical variables is evaluated in special control samples. The BABAR 
data-sample collected between November 1999 and April 2008 is used in the study of the SVT stability. 



1. Introduction 

Multi- wafer silicon (Si) tracking and vertex detec- 
tors have become an essential part of modern High 
Energy Physics experiments. Because of the short 
ionization drift distances and the sub- /im feature 
placement accuracy of silicon wafer sensors, individ- 
ual silicon wafers can provide /xm-scale position res- 
olution over areas of a few tens of square centimeters. 
In order to extend this precision over the square me- 
ter areas covered by modern detectors, the relative 
positions and orientations of the constituent silicon 
wafers must be well known (1). 

Si tracking detector construction techniques de- 
fine the wafer positions and orientations only nom- 
inally. Because Si tracking detectors are typically 
located in extremely confined spaces near the in- 



teraction region, the Si wafer positions cannot be 
measured using conventional survey techniques once 
the detector has been installed. Lab-bench measure- 
ments during construction using mechanical or op- 
tical techniques can determine wafer positions very 
accurately, but because the wafers can shift due to 
mechanical and thermal stress during and after de- 
tector installation, and because silicon charge col- 
lection effects can distort the effective position of a 
wafer from its geometric value, these measurements 
are not sufficient. Because of these effects, the wafer 
positions and orientations must be determined pri- 
marily using signals readout from the silicon detec- 
tors themselves when traversed by particles in-situ. 

This note describes the procedure developed and 
used for the BaBAR Silicon Vertex Tracker (SVT ) 
local alignment, whereby the positions and orienta- 
tions of the wafers are determined. Our procedure 
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uses track data recorded during normal BaBAR run- 
ning, filtered and prescaled to produce a fixed sam- 
ple that roughly uniformly illuminates all the wafers, 
and constrains all the local alignment degrees of free- 
dom in a statistically independent and systemati- 
cally complete way. Tracks are fit using SVT hits 
and constrained using a subset of Drift Chamber 
(DCH) and beam energy information selected to 
not impose any significant systematic bias on the 
local alignment. To avoid statistical bias, we se- 
lect an independent subset of information from each 
track. We combine track-based information with di- 
rect measurements of the relative positions and ori- 
entations of Si wafers made during detector con- 
struction, resulting in a statistically correct and sys- 
tematically robust measure of the consistency (x^) 
of a wafer's position and orientation within the de- 
tector. We use an iterative technique to determine 
the relative wafer positions that minimize the ^ 
of all wafers. The resultant local alignment is then 
validated against several possible systematic effects. 
Each of these functions are described in detail in the 
following sections. The related but simpler problem 
of determining the rigid-body position and orienta- 
tion of the SVT within the BaBAR detector {global 
alignment) is not covered in this note. 

2. The Silicon Vertex Tracker 

The BaBAR SVT was designed primarily to pro- 
vide precise reconstruction of charged particle tra- 
jectories and decay vertices near the e~^e~ interac- 
tion point of PEP-II 13), as required by BaBAR's 
diverse physics goals. Additionally, the SVT pro- 
vides the precise 6 angle measurement needed to per- 
form charged particle identification using BaBAR's 
Cherenkov detector (DIRC) and Electromagnetic 
Calorimeter (EMC) (3). 

The SVT is composed of 340 separate Si wafers, 
arranged in 5 co-axial roughly cylindrical layers, see 
Fig. 1. Each layer is composed of between 6 and 18 
modules, arranged symmetrically around the cylin- 
der (z) axis, held in place by a rigid carbon- fiber 
frame. Each module is in turn composed of between 
4 and 8 individual Si wafers, which are glued to 
supporting kevlar ribs extending in the z direction. 
See Table 1 for the geometrical parameters of each 
layer. There are six different wafer shapes, including 
a trapezoidal shape used to form the arch modules 
discussed below. The smallest wafers are 4.2 x 4.1 
cm^, and the largest (in layer 3) are 4.4 x 7.1 cm^. 



Table 1 

Geometric parameters of five SVT layers comprised of 340 
silicon wafers. The radial range for layers 4 and 5 includes 
the radial extent of the arched sections. The radius refers to 
the closest transverse radius. The length (L) and width (W) 
are along z and ij>, respectively. 
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The modules of the inner three layers are planar, 
while the modules in layers 4 and 5 are arch-shaped. 
This design reduces the amount of material and im- 
proves the point resolution for particles originating 
from the interaction region compared to a planar 
module. The modules in the inner three layers are 
tilted by 5° in azimuth (0), allowing an overlap re- 
gion between adjacent modules, see Fig. 2. This ar- 
rangement is advantageous for alignment and pro- 
vides full (j) coverage. The outer layers cannot be 
tilted because of the arch geometry. To have an over- 
lap and avoid coverage gaps in </>, the outer two lay- 
ers are divided into two sub-layers at slightly differ- 
ent radii. 

The SVT support structure is a rigid body made 
from two carbon-fiber cones, connected by a space 
frame, also made of carbon-fiber epoxy laminate. 
While in operation the SVT is mounted on the in- 
nermost magnets of the PEP-II beamline, supported 
by an assembly fixture that allows for some relative 
motion with respect to PEP-II. Because the SVT is 
mounted independently of the rest of the BaBAR de- 
tector, movement between the SVT and the rest of 
the detector can occur. During operation the SVT 
is cooled to remove the heat generated by the elec- 
tronics and is kept in a humidity controlled environ- 
ment. 

3. Parameterizing the SVT Local AUgnment 

To describe the SVT local alignment we consider 
each Si wafer as an independent rigid body located 
and oriented near its nominal construction position 
and orientation. We describe the true position of 
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Fig. 2. Transverse section of the BABAR SVT . 

a wafer by its displacement and rotation relative 
to its nominal position. The local alignment of the 
full SVT thus requires 2040 geometric parameters, 
which includes 6 redundant global degrees of free- 
dom for the position and orientation of the SVT 
within BaBAR. As described in detail in Sec. 7, we 
also model aplanar distortion of the inner 3 layers 
of wafers, adding 84 more parameters. 

We choose to define the local alignment param- 
eters in the local wafer coordinate system, a right- 
handed Cartesian system with coordinates uvw, de- 



fined uniquely for every wafer. The u axis is defined 
to lie in the nominal plane of the Si wafer pointing 
in the direction of increasing 0. The v axis is defined 
to lie in the nominal plane of the Si wafer, orthogo- 
nal to ii. The w axis is defined as the normal direc- 
tion to the nominal plane of the Si wafer, pointing 
roughly outwards from the IP. The local wafer co- 
ordinate system origin is defined as the geometric 
center of the wafer. 

The displacement component of the local align- 
ment is given by the vector {6u, Sy, dw), which de- 
scribes the true position of the Si wafer relative to its 
nominal position in that wafer's nominal local coor- 
dinate system. Similarly, we describe the orientation 
of the wafer as the vector (a^, and aw), which 
defines small right-handed rotations about the it, v 
and w axes (respectively) of the nominal local coor- 
dinate system of the given wafer, in units of radians. 
The BaBAR detector reconstruction software is writ- 
ten so that this convention of local alignment can 
be easily and efficiently applied to the reconstructed 
position of SVT hits ji). 

The readout strips on the wafers in the barrel re- 
gion of the SVT are oriented parallel to the local co- 
ordinates. The strips on opposite faces of each wafer 
are oriented orthogonally to each other, providing 
90° stereo coverage. The readout strips on the wedge 
wafers have a pitch which varies slightly along their 
length, resulting in strip directions which change 
slightly with position, but which are still roughly 
parallel. Hits reconstructed in the SVT using the 
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strips parallel to v are referred to as u hits, as that 
is the dimension they constrain. Roughly speaking, 
these hits measure the position of the traversing 
particle. Similarly, hits reconstructed using strips 
parallel to it are referred to as v hits, and they mea- 
sure the z position of the traversing particle. 

The estimated Lorentz shift in the position of u 
hits induced by BaBAR's solenoid is accounted for 
in the SVT hit reconstruction. Any difference be- 
tween this estimate and the actual Lorentz shift is 
absorbed into the Su parameter, however this can 
introduce systematic errors as described in Sec. 6. 

4. Goals and Requirements of the SVT 
Local Alignment 

The goal of the SVT local alignment is to de- 
termine the local alignment parameters with suf- 
ficient accuracy that the remaining misalignments 
contribute negligibly to the final uncertainty in the 
physics quantities extracted using the tracks recon- 
structed in the SVT . For instance, to observe CP vi- 
olation in r(45) -> B°b", the SVT must be able to 
precisely measure the roughly 250 fim average sep- 
aration between the B meson decay vertices. A full 
detector simulation study (8) showed this requires 
an average resolution of no worse than 10 //m for u 
hits and 20 /^m for v hits. To insure that the local 
alignment does not dilute these measurements, we 
require the statistical precision of the local align- 
ment contribute no more than 15% in quadrature (or 
1% net) to the vertex resolution. This implies know- 
ing 6u to roughly 1.5 /zm, 6y to roughly 3 /im, and 
aw to roughly 2 fi radians. To achieve this statistical 
precision requires a data sample equivalent to 400 
typical tracks per Si wafer. This many tracks/wafer 
is recorded in less than an hour of normal BaBAR. 
data taking. Thus meeting the required statistical 
precision is not a challenge. The real challenge of 
the local alignment procedure is to control the sys- 
tematic errors to the required level. To understand 
the issues involved in controlling the systematics we 
must first examine how the local alignment param- 
eters are constrained by data. 

A track passing through the full SVT and orig- 
inating from the interaction point (IP) will gener- 
ally generate 2 hits (1 u and 1 v) in each of 5 lay- 
ers. As a track's trajectory is well-described as a 5- 
parameter helix (^), a single track will constrain 5 
degrees of freedom in the local alignment. However, 
because tracks scatter as they pass through material, 




Fig. 3. Graphical illustration of one of global distortions that 
affect the relative position of nearby wafers only to second 
order. The elliptical effect shown has a greatly exaggerated 
scale compared to what is used for validation studies. 

the most statistically powerful local alignment con- 
straints will be on the relative positions of wafers in 
adjacent layers. Similarly, lab-bench measurements 
of relative wafer positions are useful only for nearby 
wafers, as mechanical and thermal stress uncertain- 
ties grow quickly with relative distance. 

Thus, the track-based and lab-bench measure- 
ments used by the local alignment procedure effec- 
tively only constrain the relative positions of nearby 
wafers. Many independent local constraints may of 
course be added together to completely constrain 
the local alignment, but that procedure raises the 
risk of introducing a global distortion. An example 
global distortion which correctly defines the rela- 
tive position of nearby wafers but distorts the SVT 
as a whole is shown in Fig. 3. In Table 2 we sum- 
marize the main global distortions in a system with 
cylindrical geometry, such as the BaBAR. SVT . In 
Fig. 4 we illustrate the effect of four global distor- 
tions with a natural scale of 50 fim on the position 
of individual wafers. 

Table 2 

Main systematic distortions in a system with cylindrical 
geometry and multiple layers. Distortions in r, z, and are 
considered as a function these coordinates. 
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If uncorrected, global distortions would produce 
unacceptable systematic biases in physics measure- 
ments. For instance, an uncorrected curl distortion 
would introduce an artificial charge-dependent mo- 
mentum asymmetry to reconstructed tracks, given 
the use of magnetic bending to define the charge 
and transverse momentum of a track. Likewise, a ra- 
dial scale or z-scale distortion would systematically 
change the measured distance scale of the detector, 
distorting lifetime measurements. To estimate spe- 
cific requirements on how well we must control the 
different global distortions, we consider the impli- 
cations of a subset of relevant BaBAR physics mea- 
surements. For instance, to make competitive mea- 
surements of the r lepton and B meson lifetimes, 
the absolute distance scale must be understood to 
better than 1 part in 1000. This implies controlling 
the radial scale and z-scale distortions to less than 
1/1000, or roughly 5 ^m over the size of the SVT . 
Similar arguments can be used to motivate requir- 
ing that the local alignment constrain the scale of 
all of the global distortions listed in Table 2 to bet- 
ter than 5 /im or ten times smaller than the effects 
shown in Fig. 4. The verification that our local align- 
ment procedure satisfies this requirement is given in 
Sees. 9 and 10. 

A further requirement on the the local alignment 
procedure is that it be capable of following the 
time-dependence of actual changes in the detec- 
tor. We observed some slow relative motion of the 
SVT wafers related to humidity changes, and due 
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Fig. 4. Four typical global systematic distortions out of 
nine discussed for a cylindrical system in Table 2: telescope, 
curl, bowing, and elliptical effects. Each point represents 
displacement of an SVT wafer shown with filled circles for 
layer one (•), filled triangles for layer two (A), filled squares 
for layer three (■), open triangles for layer four (A), and 
open circles for layer five (o). The typical scale of distortion 
is chosen to be 50 fim. 



to stress changes during periods of active access 
to the detector, and during changes to accelerator 
operation. The timescale for observable changes 
was about a week, which implies that the local 
alignment procedure should function on less than a 
week's accumulation of BaBAR data. 

It is also important to be able to quickly detect 
when the local alignment changes, to avoid exten- 
sive reprocessing after the initial BaBAR reconstruc- 
tion pass. To obtain feedback on possible alignment 
changes in a timely way requires that computing the 
alignment constants take no more than 24 hours. 
Because it's run frequently, the procedure must also 
be efficient in its use of computer resources. 

5. The Local Alignment Data Sample 

The data used to perform the SVT local align- 
ment are selected from those collected during nor- 
mal physics running of the -BaBAR. detector. The 
BaBAR physics trigger accepts a mix of events in- 
cluding hadronic final states of e'^e" — > qq and 
e+e^ T{AS), e+e^ l^l^ events, and cosmic 
rays which pass near the nominal IP. The IP consis- 
tency requirement of a few cm was set in the trigger 
configuration during the first three years of BaBAR 
data-taking. It was relaxed for the cosmic tracks 
with large impact parameters for studies. However, 
in our final analysis we adopt uniform approach to 
all data periods and apply the same IP requirements 
discussed below. 

Events which contain tracks useful for the local 
alignment arc identified, and the relevant tracks and 
hits are saved. These samples are passed to the min- 
imization procedure described in the next section. 
Details of the data selection are presented in the fol- 
lowing subsections, and shown graphically in Fig. 5. 

5.1. Event Pre- selection 

A pre-selection of events of eventual interest to 
the local alignment procedure is integrated into 
the BaBAR prompt reconstruction procedure that 
runs shortly after the events are recorded (lOi). All 
triggered events are first passed through a minimal 
background rejection procedure which removes most 
beam-gas interactions and scales down e'^e~ 
e+e^ interactions. A pseudo-random prescaling is 
then applied to e^e~ e+e", e+e" — > fJ,~^n~ 
and cosmic ray triggers, which results in a roughly 
uniform illumination of the detector. 



5 



track data on disk 



alignment 
constants 



ioad and fit tracks 



event seiection 







e+e- 


pairfit 




pairfit 



diagnostics 



cosmic 
join halves 



prescale in 



normal 
tracks 



hit selection 
tag overlap hits 



prescaie out 



reduced track data to disk 



Fig. 5. Diagram of the sequence for tlie event, track, and hit 
selection, including calculation of prescale factors. Sec text 
for details. 



Pre-selected events are reconstructed using the 
standard EaBAR. reconstruction program. Tracks are 
found using standard pattern recognition algorithms 
in both the SVT and DCH . Tracks foimd in the 
DCH (SVT ) are extrapolated into the SVT (DCH ) 
respectively, and hits consistent with the original fit 
are added. Tracks are fit using a Kalman filter al- 
gorithm (9) that accounts for differing hit resolu- 
tions, detector material, and magnetic field inhomo- 
geneities. 

Prompt reconstruction uses the most recent local 
alignment parameters available at the time of pro- 
cessing. If, at the end of the local alignment proce- 
dure, we observe local alignment change that might 
have affected the event selection, we repeat the pro- 
cedure using the updated local alignment param- 
eters. Selection iteration was necessary only when 
restarting the detector after a shutdown or detector 
opening. A single iteration was sufhcient to select 
an unbiased sample, because the typical changes are 
not large. The one case when several iterations were 
required was the startup of BaBAR when the start- 
ing alignment was from the optical survey. 

Pre-selected events are written to a dedicated cali- 
bration stream (file) . The calibration stream persists 
the events in the standard BaBAR mini-dst event 
format (4), which records the reconstructed tracks 
and their associated hits. In particular, this format 
records the local centroid of all SVT hits associated 
with a track. This allows measuring and applying a 
new local alignment without first having to remove 
the effect of the alignment used in prompt recon- 
struction. 



The rarest process used in the local alignment 
turns out to be cosmic rays, whose rate is inde- 
pendent of beam luminosity. Based on the cosmic 
ray trigger rate, and the number of tracks required 
to satisfy the statistical and systematic constraints, 
two days of cosmic data is sufficient to meet the 
alignment goals. We therefore define a local align- 
ment data sample based on a fixed calendar period 
of around 48 hours. 



5.2. Event Categorization and Final Selection 

A separate procedure makes a final selection of 
data useful for local alignment. This reads back the 
calibration stream, and reconstitutes the SVT and 
DCH hits using the current local alignment and cal- 
ibration. From these the full Kalman filter track fit 
is rebuilt, using the recorded hit assignments. These 
tracks are used to make a final event categorization 
and selection. 

To insure a reliable momentum measurement, 
only tracks with at least 10 DCH hits, and at least 
2 V and 3 u SVT hits (the minimum to fully con- 
strains all 5 track parameters) are selected for use in 
the local alignment. To minimize multiple scatter- 
ing effects, we also require a transverse momentum 
of at least 1 GeV/c. To cut down on background 
from secondary (material) interaction products, we 
accept only tracks whose point of closest approach 
to the BaBAR z axis is within 1.5 cm of the IP in 
the plane transverse to the z axis, and between —7 
cm and -1-9 cm of the IP along the z axis. 

Events are categorized and finally selected based 
on the multiplicity and properties of their selected 
tracks. We define three categories of events in the 
local alignment; e'^e~ — > /i^/i^ or e+e^ e+e~ 
(e^e~ l'^l~) events, cosmic ray events, and 'nor- 
mal' events. The definitions and selections of these 
categories is described below. 

Events with exactly two selected tracks are tested 
as potential e~^e~ — *■ l~^l~ events. Tracks in e"'"e~ —^ 
IJL^ l-i^ candidates are required to have associated 
EMC signals consistent with a minimum-ionizing 
particle. Tracks in e~^e~ — > e''"e^ candidates are re- 
quired to have associated EMC energy deposition 
consistent with the reconstructed track momentum. 
All e+e~ l^l~ events are required to have a total 
energy (computed from track momentum) consis- 
tent with the known combined energy of the initial 
e"'"e~ beams, and to be back-to-back in the trans- 
verse plane. 
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Candidate e+e~ — > events that pass the 

above cuts are refit using a special form of the 
BaBAR Kalman filter track fit which constrains the 
pair of tracks to have the same four-momentum as 
the initial e~^e~ system, within the independently- 
estimated errors on the beam particle momenta. If 
the of the pair fit is larger than 50, it is consid- 
ered a failed fit, and the individual tracks in these 
events are passed down to the 'normal' track selec- 
tion algorithm described below. The most common 
cause of failed pair fits is initial and/or final state 
radiation. When successful, the pair fit improves 
the track momentum resolution by more than a 
factor of 10. More importantly, the pair fit creates a 
correlated system in which information passes from 
one track through the IP to the other track. This 
allows the local alignment procedure to constrain 
the positions of wafers relative to those on opposite 
sides of the detector. We can also use e+e~ — > 
events to determine the beam momenta parameters 
simultaneously with the local alignment parame- 
ters, without using the independent beam energy 
estimate. This provides both a systematic check on 
the alignment procedure, and a precise way to mea- 
sure the beam boost. This techinque is discussed in 
Sec. 11. 

Because the BaBar track finding algorithm as- 
sumes all particle originate at or near the IP, a single 
cosmic ray passing through BaBAR is initially recon- 
structed as 2 tracks, splitting the cosmic ray trajec- 
tory through the detector roughly in half. Cosmic 
ray event candidates are selected as having two well- 
measured oppositely charged tracks which match 
in angle and position at the their point of closest 
approach to the IP. These tracks are also required 
to have associated EMC signals consistent with a 
minimum-ionizing particle. If these criteria are satis- 
fied, the hits from the upward-going track are added 
to those on the downward going track, and the com- 
bined track is refit. 

Events which fail selection as e+e~ or 
cosmic ray events are taken together with multi- 
track events as candidates for the 'normal' category. 
In this category, we select all the tracks which pass 
the usual selection cuts, plus an isolation cut. Any 
event with at least one selected track is classified as 
'normal'. 



5.3. Track Selection and Preparation 

Selected tracks from events which pass final event 
selection in any category are themselves labeled ac- 
cording to their event category. In the case of 'nor- 
mal' events, tracks are further categorized according 
to whether or not they have hits in adjacent wafers 
of the same layer in the overlap region (sec Fig. 2). 
As these overlap tracks have a very short extrapola- 
tion distance between the same-layer hits, they pro- 
vide a powerful constraint on the relative position of 
adjacent modules, and so are especially valuable in 
the local alignment procedure. Non-overlap tracks 
in 'normal' events remain categorized as 'normal'. 

To balance the impact of the largely- independent 
global distortion systematic constraints afforded by 
the different track categories, we perform a final 
track selection which roughly equalizes the statis- 
tical power of the tracks in each category for ev- 
ery individual wafer. Because the wedge wafers (sec 
Fig. 1) subtend a region of polar angle 6 where the 
rate from e+e~ — > and other physics events 

changes rapidly with 9, we further divide these into 
two roughly equal parts. 

To allow better control of the propagation of sys- 
tematic misalignment effects from the DCH into the 
SVT alignment tracks, we refit all tracks using the 
following technique. First, we split the tracks into 
two, one with all the SVT hits and one with all the 
DCH hits. Each of these associated but separate 
tracks are refit using the standard BaBar Kalman 
filter fit. The parameters and covariance matrix of 
the DCH -only track fit are sampled at the point 
where that track enters the SVT detector volume, 
and these parameters and covariance arc then used 
to constrain the SVT -only track fit. Mathematically, 
the parameter constraint is identical to the effect of 
having left the DCH hits on the track. However, by 
masking some of the parameters in the constraint, 
the information content of the DCH -only fit can be 
filtered. In particular, by masking off all but the u 
parameter (inverse curvature) of the DCH -only fit in 
the constraint, we can greatly improve the momen- 
tum resolution of the constrained SVT -only track, 
without introducing any dependence on possible sys- 
tematic distortions in the position or orientation of 
the DCH . We use the DCH -only fit uj constraint 
when fitting the e+e~ l~^l~ , cosmic ray, and over- 
lap category tracks. 

Because the sum of the local alignment param- 
eters for all wafers include the 6 global degrees of 
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freedom, the local alignment procedure could intro- 
duce a global drift. Because the uj constraint does 
not depend on the relative position or orientation of 
the SVT and DCH it cannot constrain this global 
drift. To minimize the global alignment drift, we use 
all 5 DCH -only fit parameters to constrain the fit 
of the 'normal' tracks. As these tracks have the low- 
est statistical power, this introduces only a modest 
and acceptable dependence on DCH alignment dis- 
tortions. 

5.4. Hit Selection 

From selected tracks, we select SVT hits which 
provide information useful for local alignment. First, 
hits with questionable timing or cluster shape arc 
disabled, and the tracks which held them are refit. 
Then the remaining hits are filtered to make the 
sample uniform over the detector, over several track 
categories, and over the time window in which the 
data sample was accumulated. Once the alignment 
procedure is close to convergence, a final outlier re- 
moval cut is applied. Details of the hit selection are 
shown in Fig. 5 and described below. 

For each hit on each track in each track category 
in every wafer we use a pseudo-random prescaling 
algorithm to select roughly 100 (200) hits in the 
outer two (inner three) layers, respectively, for use 
in the minimization. More hits are used for the 
inner layers to balance the larger number of wafers 
in the outer layers. The pseudo-random prescale is 
seeded on the unique event time, and so is fully re- 
peatable but effectively random. Because the data 
sample has a large number of tracks of each cate- 
gory, only a small number of hits per track are se- 
lected, reducing the number of correlated measure- 
ments used when computing the alignment x^- The 
unselected good hits are still used for the track fit. 
Cosmic events have the smallest number of tracks in 
the sample, and thus have the largest prescale fac- 
tor, corresponding to using roughly 5 hits per track. 

The hit selection is done using three passes over 
the data, interleaved with the minimization de- 
scribed in section 8. In the first pass, the hit prescale 
constants arc determined by dividing the desired 
number of hits/ wafer /category by the number ob- 
served. In the second pass, these prescale values are 
used to pseudo-randomly select hits to be used in the 
alignment x^. Selected hits are persistently tagged 
so that the same hits are used each iteration of the 
X^ minimization. In this pass we apply very loose re- 



quirements on the hit residuals, removing only the 
very worst outliers, so as to avoid biasing the align- 
ment parameters when the initial alignment is far 
from optimal. After the x^ minimization has par- 
tially converged, we repeat the hit selection proce- 
dure using the improved alignment parameters, ap- 
plying a tighter hit residual cut to suppress outlier 
hits which can distort the x^- Selected hits are again 
tagged and passed on to a final pass of x^ minimiza- 
tion. 

The hit residuals used in the alignment procedure 
are computed as the distance of closest approach be- 
tween the given hit, defined as a line in space, and 
the track trajectory, defined as a piecewisc helix in 
space, after removing the effect of that hit from the 
track fit. The residual is signed by the cross-product 
between the track direction and a nominal hit di- 
rection. The error on this unbiased residual is com- 
puted as the square-root of the quadratic sum of the 
projections of the hit error and the track covariance 
matrix onto the residual measurement. The hit error 
is estimated as a function of the hit's pulse-height 
and width, and the direction of the track. The hit 
position error functional form and parameters were 
tuned using BaBAR data. 

The hit residuals and their estimated errors are 
combined to compute a hit x^ for a particular wafer: 

hits 

xl^Y.^Vrh„ (1) 

i 

where is the residual for hit i. This x^ is relative 
to the test set of alignment parameters used when 
the track was fit. 

6. The Optical Survey Alignment 

The components of the SVT were optically sur- 
veyed to determine their relative positions at several 
stages during construction. First, individual mod- 
ules were surveyed on the lab bench during their 
construction. Then, each layer was surveyed after 
its modules were mounted on the support structure, 
starting with the innermost layer and going out. By 
studying their reproducibility, these surveys were es- 
timated to have a precision of roughly 5 /xm in the 
wafer plane and 20 /xm out of the plane. 

By averaging and combining the raw survey mea- 
surements, they were converted into a survey align- 
ment, describing the relative positions and orienta- 
tions of all the wafers in the SVT . The survey align- 
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Fig. 6. Differences in the SVT wafer positions between 
those measured by the procedure described in this paper and 
the optical survey alignment, projected in the dimensions 
illustrated in Fig. 4. Each point represents a single wafer, 
with filled circles for layer one (•), filled triangles for layer 
two (a), filled squares for layer three (■), open triangles for 
layer four (A), and open circles for layer five (o). 

ment was used as the initial condition when the lo- 
cal alignment procedure described in this paper was 
first performed. The survey alignment is also used as 
an additional constraint on relative wafer positions 
in the local alignment procedure itself, as described 
in detail below. 

A comparison of the wafer positions determined 
using the procedure discussed in this paper with the 
optical survey alignment is shown in Fig. 6. The dif- 
ferences are shown after removing the overall global 
shift and rotation between the two alignment de- 
scriptions. These figures, plus tests we made in the 
early days of BaBAR demonstrate that the survey 
alignment by itself does not meet the local alignment 
requirements defined in Sec. 4. This is understand- 
able, as mechanical stresses and other operational 



Fig. 7. Differences in the SVT wafer positions between those 
measured by the procedure described in this paper and the 
nominal geometry, projected in the dimensions illustrated in 
Fig. 4. The symbols used are defined in the caption of Fig. 6. 



effects will alter the relative wafer positions of the 
installed detector compared to lab-bench measure- 
ments. 

The survey alignment does however contain use- 
ful information. This is demonstrated by comparing 
Fig. 6 with Fig. 7, which itself compares the mea- 
sured wafer position with the SVT nominal geom- 
etry. Clearly the optical survey is a better approx- 
imation to the true local alignment than the nom- 
inal geometry. We also expect some aspects of the 
survey alignment to remain accurate even in the in- 
stalled detector. In particular, the relative positions 
of adjacent wafers in a module should be well de- 
scribed by the survey alignment, as there is little 
room for stress-induced motion between them that 
would not destroy the module. Information on the 
relative position of wafers within a module is orthog- 
onal to that provided by hit residuals, which relate 
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the relative positions of wafers in different modules. 

It is therefore desirable to add the survey alignment 
information to the local alignment procedure. How- 
ever, only the reliable parts of the survey alignment 
should be used, and the survey information must be 
appropriately combined with the hit residuals infor- 
mation. 

We add survey information to the local alignment 
procedure by constructing a survey residual for each 
wafer. The survey residual compares the wafer's po- 
sition relative to the other wafers in its same module 
as predicted by the survey alignment with the rela- 
tive position given by the test alignment, analogous 
to the hit residual definition given in Sec. 5.4. We 
assign an error to this residual based on the survey 
alignment precision, corrected for systematic effects 
which decrease the accuracy. We then construct a 
survey contribution from this residual and error, 
which is combined with the hits to form a total 
local alignment x^. 

The optical survey wafer residual is computed us- 
ing the survey measurements in the module contain- 
ing that wafer. We represent the N wafers of a given 
module {i = 1, .., A'', e.g. = 8 in layer five) spa- 
tially by a set of 9 points (j = 1, .., 9) lying in the 
wafer plane, at fixed positions in the wafer coordi- 
nate system, located roughly at the ends and mid- 
points of the sides of the rectangle defined by the 
wafer's active area. 

The wafer points in a module can be translated 
into the global BaBAR coordinate system using ei- 
ther the survey alignment, or the test alignment. We 
use the test alignment to describe the nominal global 
coordinate positions of these points . We describe 
the difference between the survey alignment and the 
test alignment transformation of these points into 
the global coordinate system as drij . 

We use the difference vectors to solve for the trans- 
lation vector R and rotation vector O which mini- 
mize the total distance between the points in global 
coordinates as described by the two alignments (test 
and survey). Together these represents the optimal 
transformation between the two alignments. The op- 
timal translation R is given by the average of the 
vector differences between the measurements 

(nxN \ /nxN \ 

where I = 1,2,3, representing the spatial coordi- 
nates. The optimal rotation O is defined implicitly 
by the equation 



3 nxN 

^^kYl i^kiinjf - {rij)k{rij)i) = 

fc=l j,i 

nxN 

where fc, I are spatial coordinates and 6ki is the Kro- 
necker delta. Solving for Q requires inverting a 3 x 3 
matrix. The components of the rotation vector Qk 
represent right-handed rotations about the respec- 
tive coordinate axis. These rotations are calculated 
with respect to the average nominal global position 
^ fij , assuming a small angle approximation (or 
equivalcntly small drij). A small number of itera- 
tions of this procedure was found to be sufficient to 
solve for the optimal transformation even when the 
rotation angles are large. This method of describing 
and solving for coordinate transforms is based on 
the formalism of rigid body rotation. It is equiva- 
lent to minimizing a constructed from drij , if we 
assign an equal error to all dimensions of all points. 

To compute the survey alignment residual, we 
consider each wafer in the module in turn, referred to 
as the wafer under consideration, with index i = I. 
For that wafer, we first calculate the transformation 
that relates the overall module position predicted by 
the survey alignment with that predicted by the test 
alignment. To avoid direct bias, we exclude the wafer 
under consideration when computing this transfor- 
mation, requiring i ^ / in Eqs. (2) and (3). To re- 
duce the impact of potential module deformation oc- 
curring after survey, we use only those wafer points 
j, i which are less than 15 cm from the center of 
the wafer under consideration when computing the 
sums. This cutoff was found to be sufficient to re- 
move systematic bias due to module deformation, 
while still providing enough points to give a statis- 
tically meaningful constraint. 

We apply the module-level transformation to 
all the points of the survey alignment, effectively 
overlaying the survey alignment on the test align- 
ment for this module's position. We then apply 
Eqs. (2) and (3) to the wafer under consideration, 
taking the transformed survey points and the test 
alignment to compute the dfij, and using only the 
points on the wafer under consideration by requir- 
ing i = I. The resulting Ai?/ and Afi/ effectively 
define a 6-dimensional survey residual, representing 
the difference between the position and orientation 
of the wafer under consideration described by the 
test alignment versus the survey alignment, relative 
to the rest of the module in the region around the 
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wafer under consideration. 

We define a 6 x 6 survey covariance matrix Vg to 
represent the intrinsic error in the survey residual. 
We use the same covariance matrix for all wafers. 
We approximate the survey covariance matrix to be 
diagonal in the local wafer coordinates, and set the 
elements according to the values found in the sur- 
vey consistency tests discussed above. Because the 
module design restricts relative wafer motion in the 
module plane, we take these values literally for the 
in-plane errors (translation in u and v, and rotation 
about w). To account for potential bowing, twisting, 
or other aplanar distortions only weakly constrained 
by the module design, we increase by a factor of 
ten the estimated errors on the remaining degrees 
of freedom (translation in w and rotations about u 
and v). 

We construct a survey from the 6-dimensional 
survey residual = {ARj , AQj) and the survey 
covariance matrix as 



X 



(4) 



The use of in the alignment procedure is discussed 
in detail in Sec. 8. 

The use of the optical survey information in 
the local alignment procedure assumes implicitly 
that the positions measured optically on the sur- 
face of the wafer correspond to the hit positions 
reconstructed in the data. In particular, the survey 
constraint could introduce a bias into the alignment 
procedure if the Lorentz shift is different for differ- 
ent wafers in a module. We have not studied this ef- 
fect in BaBAR, but we estimate it to be less than the 
estimated survey alignment errors, given the simi- 
larity of wafers in a module. Similarly, variations in 
the thickness of the wafer, which can change the ef- 
fective charge integration depth, are accommodated 
by the large out-of-plane errors we assign to Vs . 

7. Wafer Curvature 

Initial tests of the local alignment algorithm with 
BaBAR data showed a smooth but substantial vari- 
ation of residuals as a function of the local (u, v) hit 
position in some wafers, especially those in the in- 
ner layers. These variations were visible even after 
the local alignment procedure had converged. These 
effects were not seen in Monte Carlo simulation of 
the BaBAR. data or the alignment procedure. 

An example of these effects is given in Fig. 8. This 
plots a projection of the u and v hit residuals (e„ and 



ey) from a large sample of high momentum BaBAR 
tracks, as a function of the u position of the track. 
The figures show the average value of the projection, 
which is defined so as to effectively interpret the 
residuals as a local deviation in the w position of the 
wafer: 

Swu = eu/ sm{0uw) 
Swy = e„/sin(6',„u,), 

where 9uw {Ovw) is the angle between the track di- 
rection and the wafer normal in the uw (vw) plane, 
respectively. The fit to a parabola is reasonably con- 
sistent with the data, given that the errors used are 
statistical only. The points at the edges are excluded 
from the fit as they are biased by hits in the overlap- 
ping wafers in the same layer. The large uncertainty 
and fluctuations in the average u hit residuals near 
u — 1.2 cm occurs because tracks from the IP inci- 
dent at that point are nearly normal to the wafer in 
the uw plane, and so have very large error in Sw. 

An incorrect local alignment would result in a lin- 
ear dependence of Sw on u, with an offset being an 
incorrect translation and a non-zero slope being an 
incorrect rotation around the v direction. The clear 
non-linear dependence shown in Fig. 8 indicates in- 
stead an aplanar wafer distortion, not described by 
the standard six local alignment parameters. This 
geometric interpretation is supported by the fact 
that compatible effects are seen using either u or v 
residuals. We interpret Fig. 8 to say that this wafer 
is bowed in the uw plane, with a sagitta of roughly 
-40 /^m, or 15% of the wafer thickness. Bowing in 
the u direction is possible, as the support ribs con- 
strain against bowing only in the v direction. The 
large x^/NDOF (number of degrees of freedom) of 
the parabolic fits indicates that simple bowing may 
not be the only aplanar distortion present, as is dis- 
cussed further in Sec. 12. 

We see evidence for bowing in all SVT wafers. 
The observed bowing is roughly proportional to the 
u size of the wafers, with the largest effect in layer 
three. While no particular factor has been identi- 
fied which causes wafer bowing, the SVT detector 
builders agree that bowing at the observed scale is 
possible (jlll). 

The observed wafer curvature produced a sys- 
tematic bias on the transverse impact parameter as 
a function of azimuth for high-momentum tracks. 
These biases were beyond the tolerance of the local 
alignment requirements. Furthermore, because the 
wafer bowing was the same direction for all of layer 



11 



. sagitta =-45.2: 1 .8 nm 
X^/NDOF =208.7 / 43.0 




50 



u position (cm) 



^ sagitta =-38.4t 1.0 urn 
X^/NDOF =77.4 / 43.0 




u position (cm) 



Fig. 8. Average Sw projection of the u (top) and v (bottom) 
hit residuals in a layer 3 wafer as a function of the local u 
position for high-momentum tracks in the BABAR data. 



three, it caused an effective bias in the average w 
positions of roughly 30 /Ltm, well beyond the goal of 
< 5 /xm radial distortion given in Sec. 4. Thus we 
determined that the aplanar wafer distortions must 
be measured and corrected for the local alignment 
procedure to meet its requirements. 

We model the aplanar distortions as a quadratic 
dependence of the wafer w displacement 5w on the 
u position of the measurement, 



5w{u) = [u^ - ul)/2R, 



(5) 



where R is the curvature radius of the wafer, related 
to the sagitta Shyl/R = 2S/L^, L being the u half- 
width of the wafer, and Uq is a convenience param- 
eter set to Uq = L/\/3 in order to keep the w center 
of gravity of the wafer independent of the curvature 
radius R. We do not model a first-order term as that 
is redundant with the a„ alignment parameter. We 
measure R for each inner-layer wafer by fitting the 
average u and v residuals dependence on u accord- 
ing to this model, minimizing a residual-based 
to find the best R value. We do not fit for curva- 
ture in the outer layers because their smaller u size 
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Fig. 9. Results of the fits for the wafer sagitta in layer three. 
Symbols represent different modules starting from </< = in 
the order of increasing azimuth: •, I, A, A, □, o. 

makes the effect of their curvature negligible, and 
other aplanar distortions are found to dominate, as 
discussed further in Sec. 12. 

Figure 9 shows the sagitta values obtained from 
fits to the different wafers in layer three, as a func- 
tion of the wafer positions in global z. The curva- 
ture values are generally smallest near the ends of 
the module, where they are fixed to rigid hybrids, 
and are largest in the middle of the detector, where 
there is less mechanical constraint, consistent with 
expectations. 

Because the curvature parameter measurement 
depends on residuals, it is sensitive to the local align- 
ment. Similarly, the local alignment procedure de- 
pends on what value of curvature we assign to the 
wafers. This correlation forces us to fit simultane- 
ously for both the curvature parameters and the lo- 
cal alignment parameters. The organization of the 
simultaneous fit for the local alignment and wafer 
curvature is discussed in Sec. 8. We do not observe 
any time dependence to the curvature parameters, 
so those parameters are normally held fixed when 
fitting for the local alignment. 

To avoid biasing the local alignment procedure or 
BaBAR physics, we must correct for wafer curvature 
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during track reconstruction. We correct the u hits 
for curvature by displacing them in w by the 5w 
amount predicted by Eq. (5), given the measured u 
position of the hit. We correct the v hits by modehng 
them with as a three-piece piecewise hnear trajec- 
tory, where the endpoints of the three equal-length 
linear segments are chosen to lie at the Sw positions 
described by Eq. (5) , given the u coordinates of those 
endpoints. This trajectory is used when computing 
the track-hit residual, thereby naturally correcting 
for 5w. 



8. Local Alignment Minimization Procedure 

To obtain the best estimate of the true alignment 
parameters, the local alignment procedure combines 
all the available information into a total . We 
extract the optimal local alignment parameters by 
minimizing this as a function of the local align- 
ment parameters. Conceptually, this requires com- 
puting the dependance of every residual on each of 
the 2040 alignment parameters (six parameters for 
each wafer), and then minimizing the total in 
this 2040-dimensional space. However, each residual 
depends primarily on the alignment parameters of 
its hit's wafer. Additionally, minimizing such a large 
number of dimensions is computationally challeng- 
ing, raising issues of performance and accuracy. We 
therefore choose to simplify the local alignment min- 
imization procedure by dividing the total into 340 
separate wafer functions, and minimizing each in- 
dependently for that wafers alignment parameters. 
We then iterate to account for the secondary depen- 
dence of a residual on some other wafers alignment 
parameters (wafer correlation), and stop iterating 
when the alignment parameters for all wafers stabi- 
lize. 

While our iterative procedure is less direct than a 
simultaneous minimization of all parameters, we feel 
it offers numerous advantages over that technique. 
For one, standard algorithms can be used to effi- 
ciently invert the 6x6 matrices involved. Likewise, 
computing the derivatives of residuals with respect 
to a single wafer's alignment parameters is straight- 
forward and fast. Because of this, the derivatives 
can be recalculated between iterations, naturally ac- 
counting for second-order effects which are gener- 
ally ignored in a simultaneous solution. Importantly, 
our iterative procedure provides access to interme- 
diate states of the alignment, allowing us to monitor 
the convergence process directly. This gives us con- 



fidence in the final result, and allows us to test the 
sensitivity of the procedure to physical or compu- 
tational effects, as is described in Sec. 9. Addition- 
ally, our iterative procedure allows us to incorporate 
wafer curvature (Sec. 7) and beam boost (Sec. 11) 
determination during the alignment minimization, 
thus correctly handling correlations between these 
parameters and the alignment parameters. Finally, 
organizing the minimization by wafer naturally al- 
lows for a modular software design well suited for 
modern Object Oriented programming languages. 

The local alignment minimization sequence is 
shown schematically in Fig. 10. We start with the 
tracks and hits that were selected as described in 
previous sections. The tracks are fit using the cur- 
rent (test) estimate of the alignment parameters, 
from which the hit are computed for the selected 
hits. To this is added the survey alignment x^, also 
relative to the test alignment parameters, to form 
a wafer x^- We then minimize each wafer's x^ with 
respect to the change in that wafers local alignment 
parameters, holding the parameters of every other 
wafer fixed. After minimizing every wafer's x^j we 
update the alignment parameters for all 340 wafers 
by adding the computed parameter change to the 
original alignment parameters estimate. We then 
use that updated local alignment to fit the tracks 
and evaluate the survey information in the next it- 
eration, and repeat the process until the alignment 
converges. After convergence, the alignment param- 
eters are stored in the BaBAR conditions database 
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Fig. 10. Diagram of minimization sequence. Iterations over 
index TV is performed until diagnostics shows convergence. 

The wafer x^ used in the alignment minimiza- 
tion is defined as the sum of the hit and survey x^ 
defined in Sec. 5.4 and Sec. 6 respectively. We ex- 
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press that explicitly in terms of a small change in 
the wafer's alignment parameters Ap with respect 
to the test alignment used in computing the hit and 
survey residuals: 



hits 



X^^5]6f(Ap)V-ie,(Ap) 
+ ef(Ap)V7ie,(Ap). 



(6) 



We analytically minimize ^ by finding the value 
of Ap which gives a null first derivative with respect 
to all components of Ap, evaluated to first order in 
Ap. This requires computing the first derivatives of 
the hit and survey residuals with respect to Ap. The 
hit residual derivatives are calculated analytically, 
given the hit direction and the direction of the track. 
The derivative formulas are presented in Appendix 
A expressed as the Jacobian matrix: 



Jfe = aefc/c'(Ap). 



(7) 



Because the survey alignment residual is defined di- 
rectly in terms of the alignment parameters them- 
selves, its Jacobian is simply the identity 6x6 ma- 
trix. Wafers which are electronically dead or which 
have dead readout views naturally have their miss- 
ing parameters constrained by their survey align- 
ment. 

The matrix equation which results from setting 
the derivative of the wafer to zero can be inverted 
to solve for the change in the six alignment param- 
eters Ap: 



Ap 



all 

E 



■ all 

E 



(8) 



The matrix sums are computed by iterating over 
all the hit residuals and the survey residuals for 
the wafer. The new wafer alignment parameters are 
taken to be p = po + Ap, where po are the test 
alignment parameters. Upon the subsequent align- 
ment iteration, these updated parameters become 
the test parameters, and the procedure is repeated. 
The value used in convergence testing is com- 
puted directly from Eq. (6) , and thus is one iteration 
behind the parameter computation. 

A wafer is said to have converged when its total 
changes by less than a given threshold between 
iterations, typically set to 0.01 absolute. While this 
value may seem small for a that typically has a 
few hundred degrees of freedom, we found that a 
low threshold was necessary for the procedure to be 



sensitive to small global distortions, as discussed in 
Sec. 9. The entire local alignment procedure is said 
to converge when all but at most two wafers are con- 
verged. This allows for a trivial oscillation observed 
between dead or partially-dead wafers constrained 
only by survey information. The local alignment pro- 
cedure typically converges after roughly 100 itera- 
tions. 

As discussed in Sec. 7, the wafer curvature param- 
eters must be fit simultaneously with the local align- 
ment parameters. This is done in a dedicated variant 
of the local alignment procedure, where we intro- 
duce a fit for the curvature parameters between each 
normal local alignment parameter iteration, holding 
the alignment parameters fixed. The updated cur- 
vature parameters are used in the subsequent align- 
ment minimization iteration. A similar procedure is 
applied when fitting for the beam boost, as described 
in Sec. 11. 

Because both the derivative calculations and the 
X^ minimization are analytically computed, the 
minimization procedure is reasonably fast. The 
processing time for a single alignment iteration 
is limited by the time it takes to refit the tracks. 
This time is much reduced compared to the nor- 
mal BaBAR reconstruction, as the DCH informa- 
tion is applied as a single constraint instead of 40 
separate hits. A single iteration of the alignment 
minimization on a standard alignment set, without 
the curvature fit, takes roughly twenty minutes on 
a modern multi-GHz intel-processor based Linux 
computer. The entire local alignment procedure 
typically converges in roughly twelve hours. 

The local alignment minimization procedure is 
written within the i3ABAR software framework, 
using standard access to event and conditions 
data. Iterations are controlled using TCL/Tk (0), 
through either an interactive GUI or with a TCL 
script submitted in batch. Bookkeeping, diagnos- 
tics, input/output, and job management are also 
controlled through TCL/Tk. 



9. Validation of the Alignment 

We validate the performance of the SVT local 
alignment procedure through self-consistency tests 
using BaBAR data, where we compare the results 
of the aligned detector with apriori expectations. 
We also compare the performance of the aligned 
BaBAR detector with perfectly-aligned Monte Carlo 
simulation. The following tests were performed us- 
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ing data from a typical period of BABAR running, 
where the SVT aUgnment parameters were deter- 
mined according to the procedure described above, 
applying all the calibrations and corrections as nor- 
mally done when reconstructing BaBAR data. The 
validation data sample has minimal overlap with the 
data used for alignment production, due to the pre- 
scaling used in the alignment procedure. 

We use the track residuals from SVT hits to per- 
form a basic test of alignment self-consistency. The 
residuals themselves are shown in the plots on the 
left of Fig. 1 1 , with the u residuals on top and the 
V on the bottom. These show residuals from the in- 
ner three layers of the SVT , using tracks in e^e~ 
fi'^ fj,~ events selected to be within 0.2 radians of ver- 
tical in the readout-view projection. This selection 
of hits give the best resolution in the SVT , and 
so provides the most sensitivity to misalignment. 
The residual distributions are centered at zero, and 
the data and Monte Carlo distributions show very 
similar shapes. Both data and Monte Carlo distri- 
butions show some non-Gaussian tails, as expected 
from scattering tails and hit resolution variation. 
The RMS of the data and Monte Carlo distributions 
agree within a few percent. A Gaussian fit to the 
core of the distributions gives mean values consis- 
tent with zero, and sigma values of 14 (13) /Ltm for 
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Fig. 11. Single hit residuals (left) and normalized residuals 
(right) in the u (top) and v (bottom) readout view. The 
BABAR data are shown as points, Monte Carlo simulation 
as histograms. The (blue) smooth curves are the results of 
a Gaussian fit to the data. 
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Fig. 12. Difference between the impact parameter along 
the beam axis (zo); transverse to the beam axis (do), in 
azimuthal angle (</>), and in polar angle (tan(7r/2 — 6)) for 
split cosmic tracks. The BABAR data are shown as points, 
Monte Carlo simulation as histograms. The (blue) smooth 
curves are the results of a Gaussian fit to the data. 
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Fig. 13. (/>-dependence of the mean (left) and a (right) of 
the e+e~ — > fJ,'^ track cLq mismatch. The BABAR data 
are shown as points, Monte Carlo simulation as histograms. 
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Fig. 14. 9-dcpcndence of the mean (left) and cr of the (right) 
e+e~ fJ'^ /J-^ track zq mismatch. The BABAR data are 
shown as points, Monte Carlo simulation as histograms. 

data (Monte Carlo) u hits, and 18 (16) fim for data 
(Monte Carlo) v hits, respectively. The data-Monte 
Carlo width differences are consistent with the pre- 
dicted accuracy of the simulation of SVT hit reso- 
lution. 

The normalized residual (pull) distributions for 
the same tracks and hits are shown in the plots on 
the right of Fig. 1 1 , with the u residuals on top and 
the V on the bottom. Again we see good agreement 
between data and Monte Carlo. The core of the data 
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distribution is well-described as a unit-width Gaus- 
sian centered at 0. The non-Gaussian tails seen in 
both data and Monte Carlo are roughly consistent 
with known approximations in our estimate of the 
residual errors, which (for instance) do not take into 
account dead electronics channels or Moliere scat- 
tering of the tracks. 

A higher- level self-consistency test comes from fit- 
ting the incoming and outgoing branches of a cosmic 
ray as two separate tracks. Each of these tracks has 
a similar number of hits as a typical BABAR physics 
track. Both tracks are fit independently, but be- 
cause they represent the same particle, they should 
have equivalent parameters at the point whore they 
meet if the alignment is correct. In this study we 
use cosmic ray events selected as described in Sec. 5. 
Fig. 12 compares the impact parameters and angles 
of the split cosmic tracks. These tracks were fit us- 
ing their SVT hits plus a DCH curvature constraint, 
which reduces the large uncertainty in c?o and (f> com- 
ing from their correlation with curvature, which is 
poorly measured in the SVT alone due to its small 
lever arm. All tracks in these plot have a momentum 
above 2.0 GeV/c. The plots show good agreement 
between BaBAR data and Monte Carlo. Fitting the 
BaBAR data distributions to a Gaussian we extract 
single-track resolutions of 29 /xm for zq, 24 jim for 
do, 451 Iliad for cj), and 512 x 10-*^ for tan(7r/2 - 6), 
by scaling the fitted Gaussian sigma by 1/ y/2. 

Another validation test comes from comparing 
the reconstructed origin points of the two tracks 
produced in e+e~ events. Because these 

tracks are known to originate at the same point, 
the difference in their reconstructed parameters 
can be used to measure the impact parameter 
resolution, and to look for systematic biases left 
by the alignment procedure. In this study we use 
e+e" fi~^fj,~ events selected approximately as de- 
scribed in Sec. 5. To test the performance of the 
entire BaBAR tracking system, these tracks are fit 
with both SVT and DCH hits. This brings in the 
possibility that misalignments inside the DCH or 
between the DCH and SVT may affect our results. 

Figure 13 shows the ^-dependence of the e+e~ — > 
/i"*"^" tracks transverse impact parameter mismatch 
on the left, and its resolution on the right. The plot- 
ted points are the mean and the cr of a Gaussian fit to 
Edo in each (j) bin, respectively. The do mean shows 
some structure at the level of a few microns RMS, 
roughly consistent for data and Monte Carlo. We be- 
lieve this structure comes from track fit biases due 
to dead electronics in the inner layers of the SVT , 



which are partially simulated in the Monte Carlo. 

The do resolution (cr) shows a periodic variation 
due to the six-fold symmetry of the inner layers of 
the SVT (see Fig. 2), which modulates the extrap- 
olation of the hit error according to the distance 
from the innermost hit to the production point. This 
periodicity is well-reproduced in Monte Carlo. The 
Monte Carlo underestimates the do resolution by 
roughly 10%, consistent with the underestimation 
of the individual u hit residual core resolution. 

Because PEP-II produces a boosted final state, 
we cannot simply compare the e+e^ M^M" track 
longitudinal impact parameters as we did the trans- 
verse impact parameters, since the tracks are not 
back-to-back in the lab frame in the longitudinal 
projection. We can extract some information about 
the longitudinal impact parameter by constraining 
the production point to the event-average beamspot 
position, which is well measured in the transverse 
plane. However, this couples the statistical and sys- 
tematic errors of the beamspot determination with 
the alignment validation. In addition, because the 
beamspot is large (roughly 100 /xm) in the PEP- 
II bend (x) direction, the comparison has mean- 
ingful precision only for vertical tracks, where the 
beamspot constraint is limited by its measurement 
resolution of roughly 10 /xm. 

Figure 14 shows the polar-angle dependence of the 
e+e~ — > /U"'"/Lt~ longitudinal impact parameter mis- 
match on the left, and its resolution on the right. 
The plotted points are the mean and cr of a Gaussian 
fit to Azo in each 9 bin, respectively. To select verti- 
cal tracks, we use only events with track azimuthal 
angles |(/)±7r/2| <0.2. The agreement between data 
and Monte Carlo is reasonable. The observed dis- 
crepancy of the zo mean at small 9 may be due to 
systematic effects in the beamspot position determi- 
nation, which become amplified at small angles. It 
may also be related to remaining aplanar distortions 
in the outer layer wafers at large \z\, as discussed in 
section 12. The Monte Carlo underestimates the zo 
resolution by roughly 15%. This diff'erence is partly 
explainc;d by the 10% difference in intrinsic v hit 
residual core width, and by the fact that the Monte 
Carlo does not model the beamspot position mea- 
surement resolution. 

10. Validation of the Alignment Systematics 

We test the ability of the local alignment proce- 
dure to remove systematic distortions by introduc- 
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Fig. 15. Amplitude of the remaining z-expansion distortions 
as a function of iteration during the alignment procedure 
with data starting with the -50 ^tm amplitude z-expansion dis- 
tortion initial condition for data. Circles (o) illustrate proce- 
dure with survey information, while triangles (A) illustrate 
procedure with this information removed. 



ing a coherent misalignment of the SVT wafers, and 
then running the aHgnment procedure taking that 
misahgnment as the initial condition. These global 
distortions are particularly difficult to remove as the 
residuals used in the alignment procedure typically 
depend on them only to second order. We test nine 
distinct distorted initial conditions, as described in 
Table 2 and Fig. 4. We set the initial scale of these 
distortions to 50 /xm. 

In Fig. 15 we show how amplitude of the 2;-scale of 
the SVT converges back to the initial value after a 
50 /xm misalignment is applied to a standard align- 
ment set data sample. For comparison, convergence 
with and without optical survey measurements is 
shown. The survey information demonstrably pro- 
vides an important constraint on systematic distor- 
tions, which improves the convergence of the align- 
ment procedure, and reduces the systematic error of 
the final alignment. 

We fit for the amplitude of the distortion remain- 
ing as a function of iteration for each of the nine 
tests. In the data we compare wafer positions to the 
converged set of alignment parameters prior to in- 
troducing systematic distortions. In all nine cases 
we find the alignment procedure is capable of reduc- 
ing global distortions to a negligible level. We also 
perform tests with Monte Carlo where we compare 
wafer positions to the true positions known from MC 
generation. In Fig. 16 we show four representative 
initial misalignments in Fig. 4. The rate of conver- 
gence of the nine global distortions, defined as the 
decay constant of an exponential fit to the scale of 



Table 3 

Decay time (in units of iterations) for the main systematic 
distortions in the SVT local alignment procedure. The initial 
distortion is 50 fj,m in all cases. The remaining distortion is 
quoted in the table. Distortions in r, z, and (p are considered 
as a function these coordinates. 





Ar 


Az 


rA0 


vs. r 


radial 


telescope 


curl 


decay (iterations) 


5.6 


5.1 


1.3 


distortion (/im) 


0.7 


0.5 


0.1 


vs. z 


bowing 


2-scale 


twist 


decay (iterations) 


2.6 


11.2 


12.0 


distortion (/^m) 


0.6 


0.6 


0.1 


vs. <i> 


cliptical 


skew 


squeeze 


decay (iterations) 


11.8 


33.6 


32.0 


distortion (^m) 


0.9 


4.9 


4.5 



the misalignment per iteration, is given in Table 3. 

Overall, global distortions are the most weakly 
constrained deformations and it was found empir- 
ically that the order of 100 iterations were neces- 
sary to solve for these deformations. Fig. 17 shows 
the number of wafers which are not converged as 
a function of iteration. The increase after iteration 
six is due to the residual requirement applied after 
partial convergence of the procedure, as discussed 
in Sec. 8. The convergence requirement was chosen 



(a) 



(b) 



20 40 60 80 100 
Iteration 



20 40 60 
Iteration 



E 20; 



g 20 
< 



(d) 



20 40 60 
Iteration 



40 60 
Iteration 



Fig. 16. Amplitude of the remaining distortions as a function 
of iteration during the alignment procedure starting with the 
50 /im amplitude distortion initial condition : (a) telescope, 

(b) curl, (c) bowing, and (d) elliptical effects. Circles (o) 
show results with data and triangles (A) represent Monte 
Carlo. 
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Fig. 17. Number of unconvcrgcd wafers as a function of 
iteration starting with a 50 /xm amplitude elliptical expansion 
distortion initial condition relative to final alignment using 
BABAR data. 



empirically to allow convergence of the global dis- 
tortions shown in Fig. 16. 

Figure 18 plots the timc-dcpcndcncc of sonic 
global distortions of the SVT , compared to the ini- 
tial day-one alignment. The plot covers roughly 50 
time periods when the internal structure of the SVT 
was suspc;ctcd to be changing due to mechanical 
stress during detector access or humidity changes. 
The 2/-axis plots the amplitude of the change in a 
particular distortion, obtained using a method sim- 
ilar to that shown in Fig. 4. The number of days 
since the initial alignment is shown on the a;-axis. 
The first two points at negative time compare the 
initial alignment to the ideal geometry (day —300) 
and the survey geometry (day —200), as shown in 
Figs. 7 and 6 respectively. The large bowing in the 
outer layers around day 700 is due to an accidental 
humidity increase, which caused the carbon-fiber 
support structure to expand. 

Our study of potential distortions of the SVT 
wafer positions in the alignment procedure places 
limits on systematic uncertainties in physics mea- 
surements, such as particle lifetime or the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix parameter 
measurements through time-evolution studies of 
the B meson decays. There are two sources of sys- 
tematics: any misalignment due to time- variations, 
such as those shown in Fig. 18, and due to imper- 
fections in the alignment procedure with residual 
misalignments remaining, both statistical and sys- 
tematic. We minimize the former by having about 
50 independent time periods which follow major 
changes in the detector. The systematic distortions 
are controlled in the validation plots, such as those 
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Fig. 18. Timc-dcpcndcncc of the global distortions of the 
SVT when compared to day-one. Each point represents an 
alignment set and the day since the first alignment set is 
shown on the x-axis, while days —300 and —200 correspond 
to ideal and survey geometries. Amplitudes of four distor- 
tions are shown: (a) telescope, (b) curl, (c) bowing, and (d) 
elliptical effects. In (a) and (b) triangles (A) represent aver- 
age effects for all wafers, while in (c) and (d) the inner three 
layers and the outer two layers are shown separately with 
the filled (•) and open (o) circles, respectively. 

shown in Figs. 12 and 13. We have validated that 

all major systematic distortions shown in Table 2 
would be visible in the above validation plots. 

To facilitate studying the impact of potential re- 
maining misalignment on physics analysis, we cre- 
ate special alignment sets with intentional misalign- 
ments, which describe remaining misalignments pos- 
sible by either of the mechanisms above. The impact 
of remaining misalignment on a physics analysis is 
evaluated by reconstructing the tracks in simulated 
events using a misaligned parameter set, followed by 
the normal analysis chain. The difference between 
that result and the one produced from analyzing the 
same events with perfectly aligned track reconstruc- 
tion is taken as the misalignment systematic error. 
The misalignment error is rarely the dominant sys- 
tematic error of a BaBAR analysis. 

11. Fit for the e+e" Beam Energies 

In the joint fit of the e+e^ the same 

four-momentum for the initial state (from the beam 
energy monitoring) and the final state is assumed. 
However, while beam angle measurements were 
found to be precise in the PEP-II beam monitor- 
ing, measurements of the beam energies were not 
stable to better than a couple of MeV. In addition. 
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Fig. 19. Distribution of the dimuon /i^^^ invariant mass 
(left) and total measured momentum of the two muons 
(right) in data (points with error bars) and MC (histogram). 

initial and final state radiation can systematically 
change the total energy and momentum of the 
e+e~ — * Z+/~pair compared to the original e+e~. 
Assuming the wrong momentum constraint in the 
pair fit could introduce a telescope distortion, as 
described in Table 2. This effect was indeed ob- 
served in early local alignment validation. We use 
the cosmic tracks to constrain this effect. Because 
these cross the entire detector, they would detect a 
telescope distortion as a kink in the track polar an- 
gle going from one side to the other. Because of the 
cosmic track constraint, we can use the SVT local 
alignment procedure to determine simultaneously 
the alignment and the boost of the two muons. 

In Fig. 19 we show the distributions of the dimuon 
invariant mass and total measured momentum of 
the two muons (boost). In the alignment procedure, 
we fit the average boost and use it in the 

four-momentum pair constraint instead of the initial 
state boost. We then iterate, allowing a simultane- 
ous extraction of the boost value and the geomet- 
rical constraints on the telescope global distortion. 
This simultaneous fit makes the local alignment al- 
gorithm more stable, and allows for monitoring the 
e+e~ beam boost. In principle, we can also use the 
e~^e~ events to measure the total e~^e~ in- 

variant mass. However, because of the high momen- 
tum, and therefore poor relative momentum resolu- 
tion, of the tracks in e+e^ events, the reso- 
lution is poor. A better estimate of the e+e~ invari- 
ant mass can be made using T(45) — > BB decays, 
when one B meson decays fully hadronically. 

12. Residual Aplanar Distortions 

As described in Sec. 7, we see clear cvidcincc; for 
substantial aplanar distortions of wafers in the in- 
ner layers of the SVT . The dominant effect can be 
characterized as a bowing of the wafers along their u 



coordinate, which has limited external support. We 

also studied more general forms of aplanar distor- 
tions using the same techniques discussed in Sec. 7, 
where we interpret the average residuals as due to 
a local Sw distortion. From these studies we discov- 
ered that more general forms of aplanar distortion 
are present in the BaBAR SVT . 

An example is shown in Fig. 20. This plots the 
average Sw of layer five v hit residuals for high- 
momentum tracks in BaBAR data as a function of 
their position. The distortions arc substantial, but 
they cannot be described by simple u bowing of in- 
dividual wafers. Some patterns are evident, such as 
the V bowing of the wafers centered around z = — 10 
cm. A similar pattern is seen at the forward end of 
layer five, and at both ends of layer four, but not in 
similar locations in layers 1^3. This suggests that 
these distortions may be related to the bend at the 
end of the arch modules (see Fig. 1). Fig. 20 also 
shows irregular distortions in many layer five wafers 
near their edge at approximately z = 5 cm. This 
is where the forward and backward module halves, 
which were constructed separately, join. Similar dis- 
tortions occur where the layer four module halves 
join. No such distortions are seen in the middle of the 
inner layer modules, which were built in one piece. 

The aplanar distortions in the outer layers are 
thought to be responsible for some of the remain- 
ing irregularities seen in the validation plots, for in- 
stance the variation of e+e~ miss distance 
resolution with (p and 9, as shown in Figs. 13 and 
14. Unfortunately, the irregularity of these distor- 
tions makes them difficult to correct, and we do not 
attempt to do so in the BaBAR local alignment pro- 
cedure. A full Sw map of a vertex detector would in 
principle be possible given a large and diverse col- 
lection of data. In particular, tracks originating at 
many different positions would be necessary to avoid 
the lack of constraint on Sw when the projected track 
incident angle is normal to the wafer. 

13. Conclusions 

We have described the procedure used to deter- 
mine BaBAR SVT local alignment. We have shown 
that this procedure satisfies the requirements placed 
on the SVT performance by the BaBAR physics 
goals. We have demonstrated that this procedure 
is robust against global distortions that could oth- 
erwise introduce unacceptable systematic biases in 
BaBAR tracking and physics data. 
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Fig. 20. Average <5to projection of the v hit residuals in the backward half of layer five as a function of the global if> and z 
position for high-momentum tracks in the BABAR data. The overlaid lines show the approximate wafer edges. The forward 
half of the module near 4> = —0.25 radians shows no data due to its failed v view readout electronics. 
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Appendix A. Derivative calculations for the 
SVT Local Alignment 

To minimize the used in the SVT local align- 
ment procedure, we need the first derivatives of the 
track-hit residuals with respect to the six degrees of 



freedom of the wafer alignment defined in Sec. 3, be- 
ing specifically the three translations along the local 
wafer coordinate axes {u, v, and w) and the three 
small rotations about those axes. In BaBAR, residu- 
als are defined as the distance in space between the 
track and the hit trajectory at their point of closest 
approach (FOG A), signed by the cross-product of 
the track direction (t) and the hit trajectory direc- 
tion (h). The hit trajectory is defined as a fine seg- 
ment in the wafer plane, with a direction given by 
the strips used in this hit (generally u or v ). Thus 
a barrel-module u hit (which constrains the u posi- 
tion of a track) has hit trajectory direction h — v, 
and a v hit has hit trajectory direction h — u. 

The derivatives of a residual e with respect to 
wafer translations d and rotations a can be ex- 
pressed as: 



de_ 

da. 



(A.l) 
(A.2) 



where D = t x h/\t x h\, H is the position of the 
hit at FOGA, relative to the geometric center of 
the wafer, and i € {u,v,w}. These derivatives are 
coded in the BaBAR local alignment procedure using 
the CLHEF (jl2r ) class library geometry methods. 
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The exact derivative calculations and implementa- 
tion were tested by comparison with numerically- 
computed derivatives. 
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